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SECTION  I 


INTRODUCTION 

^Stable  high-density  plasma  discharges  have  many  applica¬ 
tions  in  plasma  physics.  An  especially  significant  example  is 
the  use  of  these  discharges  in  high-power  gas  lasers.  The  work 
described  here  concerns  investigation  of  the  external  ionization 
e-beam  sustainm  method  of  producing  a  high-density  discharge, 
together  with  the  advantages  and  instabilities  peculiar  to  that 
system.  _ 

The  high-pressure  electron  beam  stabilized  (KBS)  discharge 
exists  stably  in  two  regimes,  distinguished  by  their  modes  of 
electron  production  and  loss.  First,  if  secondary  ionization  is 
much  less  than  EB  ionization  (i.e.,  if  S  •  AiVenc  where  S  is  the 
ionization  source  function,  Aj  is  the  first  Townsend  coefficient, 
ne  is  electron  number  density  and  V0  is  electron  drift  velocity), 
we  have  the  conventional  EBS  discharge,  in  winch  electron  produc¬ 
tion  is  balanced  by  electron-neutral  attachment,  by  electron¬ 
positive  ion  recombin.it  ion,  or  in  general  by  a  combination  of 
both.  Second,  if  S  AfVene  we  have  the  so-called  avalanche 
mode;  the  EB  acts  as  a  trigger,  and  the  number  density  of  elec¬ 
trons  increases  until  the  number  lost  by  recombination  equals  the 
number  gained  by  Townsend  multiplication.  We  note  that  in  gases 
exhibiting  electron  attachment,  there  is  a  third  mode  of  discharge 
operation  in  which  negative  and  positive  ion  current  is  much  great¬ 
er  than  electron  current.  In  the  work  described  here,  the  dis¬ 
charge  operates  in  the  first  mode,  and  recombination  processes 
dominate  the  electron  losses. 

1.  ATTACHMENT  INSTABILITIES 

Attachment  instabilities,  discovered  and  investigated  under 
this  program, 1«2  occur  when  the  differential  gas  conductivity  be¬ 
comes  negative  due  to  enhanced  electron  attachment  by  neutral 
molecules  at  higher  electric  fields,  and  are  significant  at  high 
fields  and  low  pressures  m  many  gases.  Thermal  and  acoustic  m- 
stabilities^  occur  generally  for  discharges,  providing  limits  on 
the  total  energy  that  can  be  introduced  into  the  discharge.  As 
pointed  out  by  Nighan  and  Wiegand,4  vibrational  instabilities 
exist  when  significant  energy  is  stored  in  molecular  vibrational 
modes;  increase  in  gas  temperature  accelerates  the  collisional 
transformation  of  this  energy  to  translational,  resulting  in 
runaway  heating  of  the  discharge  medium. 

1.  Douglas-Hami  1  ton,  D.H.,  Mam,  S.A.,  J.  Appl.  Rhys.  4J>,  4406  (  1  974) 

2.  AFOSR  Final  Report,  Contract  F44620-70-C-0023  (1974). 

3.  Jacob,  J.H.  and  Mam,  S.A.,  Appl.  Phys.  Lett.  26,  53  (1  975). 

4.  Nighan,  W.L.  and  Wieqand,  W.J.,  Appl.  Phys.  Lett.  25,  633  (1974). 
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2. 


STREAMERS 


Both  attachment  and  therma 1 -acoustic  instabilities  are 
large-scale  positive  column  instabilities,  and  only  the  initial 
phases  of  their  growth  can  be  described  by  linear  analysis.  The 
secondary  phases  of  all  these  instabilities  appear  to  be  the 
formation  of  discharge  streamers*,  and  we  are  now  beginning  to 
obtain  detailed  information  about  this  nonlinear  phase  of  the 
discharge  instability. 

Streamers  of  this  type  in  EBS  discharges  were  first  dis¬ 
covered  at  AERL  under  this  program.  Since  that  time,  streamers 
have  been  found  to  occur  quite  commonly  in  such  discharges,  and, 
indeed,  are  responsible  for  arcing  in  many  high  power  gas  dis¬ 
charge  lasers.  Examples  are  shown  in  Figures  1  and  2.  Accord¬ 
ingly,  in  recent  work  under  this  program,  we  have  concentrated 
on  study  of  the  streamer  phenomenon  and  have  succeeded  both  in 
verifying  the  essential  elements  of  our  original  model  of  a 
streamer,  and  in  obtaining  new,  more  detailed  data  about  streamer 
propagation . 

In  our  original  model,  shown  in  Figure  3,  a  streamer  was 
described  as  a  cylinder  of  hot  and  consequently  conductive  gas 
(both  because  of  its  lower  density,  and  because  of  the  thermal 
ionization  of  the  gas  at  the  equilibrium  streamer  temperature 
near  6000°K).  This  conducting  cylinder  protrudes  into  the  dis¬ 
charge;  near  its  head  there  will  be  both  a  strong  field  and  a 
high  current  density  (just  as  would  occur  near  the  tip  of  a  wire 
protruding  from  the  cathode) .  Both  field  and  current  act  to 
raise  the  conductivity  of  gas  beyond  the  tip  -  the  field  by  ioni¬ 
zation  and  the  current  by  heating.  The  streamer  grows  across  the 
discharge  at  a  velocity  determined  by  the  time  required  to  heat 
the  gas  up  to  conducting  temperature.  When  it  reaches  the  appro¬ 
priate  electrode,  or  joins  with  a  streamer  propagating  from  the 
opposite  electrode,  a  low- impedance  circuit  is  formed  and  an  arc 
interrupts  the  discharge. 


•Non  self-sustaining  discharge  streamers  must  be  distinguished 
from  the  electron-avalanche  type  of  streamer  appearing  in  spark 
breakdown.  In  the  first  case  the  streamer  closely  resembles  an 
incomplete  arc  and  photioni zat ion  processes  may  be  neglected, 
where  in  the  second  case  phot  10m zation  generates  the  seed  elec¬ 
trons  in  the  avalanche.  In  the  present  work  we  discuss  only  the 
non  self-sustaining  discharge  streamer. 
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ANODE 

CATHODE  RODS 


Figure  1  Pul  :.•••!  Diseh-ii  :•  Streamers.  Flection  density 

,  electric  field  5000  v/cm, 
p re:- sure  1  atm.  > i :■  mixture  He:N2?C02»  3:2:1. 
Electron  : ■  \r  enters  through  cathode,  at  bottom. 
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SECTION  II 


EXPERIMENT,. L  RESULTS 


The  first  time-resolved  photographs  of  streamers  were  ob- 
Mim\!  the  previous  year.  Using  a  high-speed  (44,000  frames/ 
second  framin;  camera,  we  were  able  to  record  successive  stages 
1 r  the  growth  f  a  single  streamer,  which  provided  us  with  much 
hotter  data  about  the  structure  and  the  propagation  of  the 
streamer.  The  photo  sequences  showed  that  streamers  do  indeed 
have  the  character ist ic  shape  of  a  cylinder  with  a  ball  on  the 
end.,  is  had  been  :  redicted  Ly  our  theoretical  model,  and  that 
within  experimental  accuracy,  the  growth  of  a  streamer  is  an 
exponential  process.  Moreover,  the  exponential  growth  rate 
deduced  from  the  photo  sequences  correlates  well  with  the  arcing 
time  .ten  in  the  discharge  current  traces. 

Data  were  taken  in  nitrogen  at.  atmospheric  j  ressure  over 
i  wide  range  of  discharge  currents  and  electric  fields.  When 
*  he  lata  were  rrouped  into  slow,  intermediate,  and  fast  growth 
siscs,  and  displayed  parametrically  in  the  cur  rent -power  plane, 
i  clear  dependence  was  apparent.  As  shown  in  Figure  4,  a  straight 
line  can  bo  drawn  in  this  plane  to  separate  the  slow  growth 
ases  from  the  intermediate  and  fast  growth  cases.  This  "line 
f  constant  rrowth  rate”  corresponds  to 

P  -  3.8  j  -  1  (1) 

where  P  is  in  kW  cm-*  and  j  is  in  A  cm^.  As  can  bo  seen  from 
Figure  4,  this  line  is  quite  precisely  defined  by  our  data  over 
a  considerable  parametric  ranae. 

In  the  P- 1  plane,  a  straight  line  through  the  origin  would 
represent  a  constant  value  of  electric  field.  The  observed 
instability  line  is  almost,  but  not  quite  of  this  form.  The 
difference  is  apparent  when  the  data  are  plotted  against  electric 
field.  This  is  done  in  Figure  5,  where  the  same  data  points, 
and  the  same  separation  line  are  shown  in  the  P-E  plane.  Here 
the  line  is  a  hyperbola.  At  large  powers,  the  line  is  nearly 
vertical,  implying  that  *  is  there  dependent  primarily  upon  F. 

Put  at  lower  powers  and  hiqher  electric  fields,  the  line  is 
nearly  horizontal,  indicating  that  is  mainly  P  dependent  in 
that  region. 

In  the  early  work  just  described,  the  framing  camera  was 
focused  on  a  small  strip  of  the  discharge  extending  from  cathode 
to  anode,  and  streamers  were  made  to  form  there  by  attaching  a 


(  ) 


Figure  5  Dependence  of  Streamer  Growth  Bate  Upon  Discharge 
Power  Density  and  Electric  Field.  y  is  in  units 
of  (100 


small  loop  of  tunqsten  wire  to  the  cathode.  This  forced  the 
streamers  to  start  with  a  size  of  about  1  mm  and  to  always  start 
at  the  same  place,  which  made  the  data  reproducible,  permitting 
us  to  determine  Eq .  (1),  but  which  prevented  us  from  learning 

anything  about  streamer  initiation,  or  comparing  streamer  growth 
at  different  points  within  the  discharge. 

During  this  past  year,  we  improved  our  experimental  appa¬ 
ratus  and  technique  to  the  p>oint  where  we  could  obtain  clear  and 
detailed  framing  camera  photographs  of  the  entire  discharge. 

This  enabled  us  to  dispense  with  the  tungsten  wire  and  study 
spontaneous  growth  of  streamers  from  a  smooth,  rounded  cathode 
plate.  The  rebuilt  apparatus  is  shown  schematically  in  Figure  fe. 
The  actual  device  is  shown  in  Figures  7a  through  7d.  It  includes 
a  completely  new  discharge  cell  which  is  cleaner,  more  convenient, 
and  affords  a  much  improved  view  of  the  discharge  region.  The 
cell  contains  a  reflecting  prism  which  gives  a  side  view  of  the 
discharge  through  a  window  in  the  end  plate  of  the  cell.  We  also 
added  to  the  cell  a  pair  of  coils  and  associated  pjulsing  circuitry 
that  can  be  used  to  impose  an  external  maqnetic  field  transverse 
to  the  discharqe  current.  This  field  exerts  a  sideways  force  on 
the  discharge,  but  for  fields  less  than  200  G,  this  force  should 
not  have  much  effect  on  a  uniform  discharqe  in  our  apparatus. 
However,  in  a  streamer  the  current  density  is  orders  of  magnitude 
above  the  background  level  and  our  calculations  suggested  that  the 
proportionally  larger  force  miqht  be  significant.  A  sideways 
force  would  tend  to  push  the  streamer  through  the  surronding  cold 
gas,  which  should  help)  slow  down  its  growth.  Addition  of  the  coils 
enabled  us  to  determine  whether  fields  compatible  with  e-beam 
piropiagat  ion  could  have  an  observable  effect  on  streamer  growth. 

In  this  year's  experiments,  the  framing  camera  was  run 
in  half-frame  mode  (22,800  f rames / second )  to  obtain  a  longer 
expiosure  and  a  larger  f  leld-of-view.  A  small  resistance  (5  .. ) 
was  added  in  series  with  the  discharge  to  reduce  the  brightness 
of  the  final  arcs  which  had  caused  back-fogqing  of  many  of  the 
framing  camera  sequences  taken  the  previous  year.  These  changes 
enabled  us  to  obtain  much  better  photoqrap>hs  of  streamers. 

Examples  of  the  new  data  are  shown  in  Figures  8  through  13 
which  show  streamers  in  nitroqen  discharges  at  various  pressures, 
electric  and  magnetic  fields,  and  power  loading.  Each  frame  of 
each  sequence  shows  the  entire  width  of  the  discharge  and  the 
f leld-of-view  extends  from  the  cathode  to  about  two-thirds  of  the 
way  to  the  anode.  The  cathode  is  at  the  bottom  of  each  frame. 

(It  can  be  seen  in  some  after  the  arc.)  The  view  is  exactly 
across  the  face  of  the  cathode.  The  framing  rate  is  one  frame 
every  4  4  „s  and  the  exposure  of  each  is  18  ,.s.  The  e-beam 
entered  through  the  anode  (a  plate  with  a  hexaqonal  pattern  of 
holes)  and  the  maqnetic  field,  when  present,  is  along  the  line- 
of-sight.  The  width  of  the  cathode  is  5  cm  and  the  anode  to 
cathode  spacing  is  1  in. 
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Figure  6  Schematic  Diagram  of  BBS  Facility 
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Figure  10  Streamer  Growth 
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For  each  such  shot,  oscilloscope  traces  of  the  discharge 
voltage  and  the  discharge,  e-beam,  and  magnet  coil  currents 
were  recorded,  as  well  as  the  composition  and  pressure  of  the 
fill  gas,  the  e-beam  current  setting,  and  the  charging  voltages 
of  the  sustainer,  e-beam,  and  magnet  coil  banks.  We  thus  have 
a  quite  detailed  description  of  each  shot,  as  shown  in  Table  1 
which  lists  the  data  for  the  shot  shown  in  the  photo  sequence 
in  F igure  1 3 . 


Figure  8  shows  a  typical  discharge  which  produced  an  arc 
155  s  after  the  beginning  of  the  pulse.  Many  streamers  can 
be  seen  and  the  last  frame  clearly  shows  which  one  produced  the 
arc.  Figure  9  shows  a  similar  shot  at  a  slightly  lower  voltage. 
Mere  the  time  to  arc  was  longer  (265  ..s),  so  the  streamer  growth 
can  be  seen  in  more  detail.  A  whole  row  of  streamers  is  evident 
growing  from  one  edge  of  the  cathode.  It  is  interesting  to  note 
the  bright  spots  at  the  edges  of  the  cathode.  These  are  grooves 
in  the  side  of  the  cathode  plate  (see  Fiqure  14)  which  appear  to 
be  acting  as  hollow  cathodes. 


Figure  10  shows  a  discharge  under  conditions  identical  to 
Figure  9  except  for  the  addition  of  an  external  magnetic  field. 
The  B-field  had  no  clear  effect  on  streamer  growth,  and  no 
sideways  motion  of  the  streamer  is  apparent.  But  note  that  this 
discharge  arced  in  the  center,  while  in  Figure  9  the  arc  is  at 
the  edge.  This  is  a  surprising,  but  reproducible  result.  On 
.ill  six  runs  we  took  with  an  imposed  B- field,  the  arc  was  in 
the  center.  In  ten  runs  in  nitrogen  without  a  B-field,  the  dis¬ 
charge  usually  (seven  times)  arced  at  the  edge.  It  should  be 
noted  that  an  arc  which  appears  to  be  in  the  center  might  be  on 
the  front  or  back  edge;  but  in  Figure  11,  which  shows  the 
cathode  after  completion  of  these  experiments,  one  can  see  burn 
marks  where  the  arcs  occurred,  and  some  are  in  the  center. 
Indeed,  there  are  two  groups  of  burn  marks,  one  set  along  the 
edges  and  another  set  clumped  in  the  center,  with  only  a  few 
intermediate  cases.  Evidently,  this  distribution  is  related  to 
the  streamer  initiation  process  which  we  plan  to  explore  further 
in  subsequent  work. 


Figure  11  shows  another  new  observation:  a  filmed  record 
of  a  post  arc.  In  this  shot,  the  e-beam  shut  off  before  the 
streamer  had  reached  the  anode  and  the  discharge  current  dropped 
slmost  to  zero.  But  the  streamer  continued  to  grow  and  produced 
an  arc  roughly  100  „s  later.  Figure  12  shows  an  even  more  pro¬ 
nounced  example.  This  arc  did  not  occur  until  about  a  milli¬ 
second  after  the  discharge,  but,  as  can  be  seen,  it  was  still 
caused  by  formation  of  a  streamer.  These  films  are  our  first 
direct  evidence  that  streamers  are  responsible  for  post  arcs. 

Figure  13  shows  a  discharge  in  which  several  streamers 
grew  to  significant  lengths  before  the  strongest  one  produced 


TABLE  1 


DATA  FROM  RUN  9,  OCT.  19,  1978  (Film  16) 


Pressure  .ind  gas:  500  torr  N2 
Sustainer  bank  voltage:  8  kV  , 

e-beam  current:  0.28  A,  0.67  raA/cm*' 

e-beam  voltage:  90  kV  unattenuated,  approx.  75  kV  after  foil 

Coil  bank  voltage:  1.0  kV 

Magnetic  field:  lt>b  gauss 

Apparent  location  of  arc:  center 

e-beam  pulse  duration:  220  v.s 

Time  to  arc:  227  i.  s 

Streamer  exponentiation  Time:  86  s 
e-foldings  to  arc:  2.56,  -2.64 


I  n  1 1  1  a  1 

Final 

Sustainer  voltage  (kV) 

7.66 

5.66 

Electric  field  (kV/cm) 

3.01 

2.23 

Sustainer  current  (A) 

36.0 

32.0 

Current  density  (A  cm?) 

1.44 

1 .28 

Power  leposition  (kW  cm  -  ! 

4.34 

2 .85 

Dischar  :e  resistivity  (k..-cm) 

2.09 

1  .  74 

E  N  ( kv  cm- a to 

4 .58 

3.  39 

Energy  deposition  (J  cm  J) 

.81 

Comments:  At  least  seven  sizable  streamers  arc  visible.  In  the 

early  frames  those  near  the  edge  are  larger,  but  are 
later  overtaken  by  a  brighter  streamer  near  the  center. 
The  last  frames  also  show  a  second  streamer  near  the 
center  which  is  relatively  bright  for  its  size  and 
appears  to  be  growing  rapidly. 


Showing 


A  View  of  the  Cathode  After  the  Experiments 
Burn  Marks 


Figure  14 


an  arc.  Here,  as  in  some  of  the  other  films,  one  can  see  stream¬ 
ers  bending  toward  the  center  of  the  discharge.  It  is  not  clear 
if  this  is  due  to  bending  of  the  electric  field  or  to  a  gradient 
m  the  background  conductivity,  current  density,  and  gas  tem¬ 
perature.  (The  e-beam  is  nonuniform,  being  peaked  at  the  center 
due  to  scattering  on  passage  through  the  foil.) 

Films  were  also  taken  of  shots  in  3:2:1  laser  gas,  but 
although  these  discharges  did  arc,  streamers  are  not  seen  in 
the  photographs.  As  can  be  seen  in  Figure  15,  glowing  spots  do 
appear  on  the  cathode  early  in  the  discharge,  and  a  conical 
structure,  presumably  the  base  of  a  streamer,  is  visible  on 
the  cathode  }ust  before  the  arc,  but  the  streamers  are  apparently 
too  cool  to  radiate  sufficiently  in  the  wavelength  range  to 
which  the  film  is  sensitive.  This  is  not  too  surprising,  since 
a  lower  temperature  is  sufficient  to  ionize  CO2 .  Streamers 
have  been  seen  in  3:2:1  gas  in  other  discharge  cavities,  using 
integrated  exposure  techniques  and  exposing  the  gas  resistant 
discharges.  In  the  future,  we  hope  to  use  shadowgraph  techniques 
to  observe  streamers  of  this  sort. 

Experiments  were  also  done  in  helium,  but  the  streamer 
mode  of  arcing  was  not  seen.  The  conventional  breakdown  thresh¬ 
old  is  so  low  in  helium  that  a  voltage  that  did  not  make  an  arc 
before  the  e-beam  pulse  was  also  insufficient  to  produce  streamer 
growth  during  the  few  hundred  microseconds  of  a  discharge. 

Altogether,  this  year's  work  has  produced  photograp'hic 
records  of  25  such  runs  as  well  as  a  much  larger  number  of 
osc 1 1 lograph  records  of  runs  in  which  the  frame  camera  was  not 
used.  A  listing  of  the  runs  which  do  have  framing  camera 
records  is  given  in  Table  2. 

The  new  data  on  nitrogen,  which  now  includes  discharges  at 
low  pressures,  are  consistent  with  our  earlier  results.  We  find 
that  streamer  growth  is  exponential,  except  in  the  early  stages 
when  the  streamer  is  less  than  a  few  millimeters  long.  In  tha^ 
initial  stage,  which  was  not  seen  last  year  because  of  the  use 
of  the  loop  of  tunqsten  wire,  the  growth  appears  to  be  faster 
than  the  exponentiation  seen  later.  We  also  find  less  correla¬ 
tion  now  between  the  time  to  arc  and  the  observed  growth  rate 
of  a  developed  streamer.  There  is  some  variability  in  the  time 
required  to  initiate  a  streamer,  an  observation  which  gives  cause 
for  optimism  that  a  proper  choice  of  cathode  can  reduce  the  fre¬ 
quency  for  arcing  in  such  discharges.  In  the  later  stages  of 
growth,  we  observe  an  exponentiation  time  with  parametric  scaling 
like  that  seen  before,  provided  that  the  electric  field  and  the 
power  density  are  both  scaled  to  the  pressure  of  the  gas. 
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SECTION  III 


THEORETICAL  ANALYSIS  AND  COMPARISON  WITH  DATA 


1 .  STREAMER  MODEL 

Our  original  model  of  a  streamer,  shown  in  Figure  3  above, 
included  a  core  cylinder,  essentially  an  arc,  of  length  c  and 
radius  b,  a  cap  sphere  or  halo  of  radius  a,  where  the  conductivity 
is  high  due  to  ionization  by  the  field  and  current  concentration 
around  the  tip  of  the  streamer  column,  and  finally  the  background 
discharge  of  cold  gas  with  uniform  conductivity  produced  by  the 
e-beam.  This  model  was  postulated  before  we  had  any  detailed 
data  about  streamer  structure,  but  it  has  been  well  borne  out  by 
our  subsequent  experimental  work.  As  is  evident  in  the  framing 
camera  photographs  in  the  preceding  section,  streamers  quite  con¬ 
sistently  display  the  structure  which  we  had  predicted. 

In  previous  work  under  this  contract,  the  model  was  quanti¬ 
fied  in  considerable  detail  and  used  to  predict  such  parameters 
as  the  temperature  and  the  diameter  of  the  streamer  column,  the 
size  of  the  luminous  halo,  and  the  rate  of  growth  of  the  streamer 
under  various  discharge  conditions.  The  predicted  streamer 
growth  rate  is  in  partial  agreement  with  our  observations.  The 
streamer  velocity  deduced  from  the  model  is  certainly  of  the  right 
order,  which  gives  us  confidence  that  the  basic  mechanism,  resis¬ 
tive  heating  of  the  gas  near  the  streamer  tip,  is  correct.  As 
reviewed  below,  the  model  predicts  that  the  streamer  growth  rate 
should  be  primarily  dependent  on  the  power  density  in  the  back¬ 
ground  discharge.  This  is  found  to  be  true  in  some,  but  not  all 
of  the  parametric  regimes  which  have  been  studies  in  our  e-beam 
discharge. 

To  improve  the  accuracy  of  the  model  we  bel leve  it  neces¬ 
sary  to  describe  in  more  detail  the  heating  processes  around  the 
streamer  tip.  The  tip  is  a  complicated  structure  where  the  gas 
and  electron  densities  and  temperatures  and  the  electric  field 
strength  and  current  density  are  all  rapidly  varying.  Beyond 
the  tip  the  conductivity  is  due  to  e-beam  ionization  of  the  gas. 
Within  the  halo  we  believe  the  conductivity  is  due  to  ionization 
by  electron  multiplication,  because  the  field  is  near  the  criti¬ 
cal  field,  and  to  thermal  ionization  since  the  gas  is  heated  by 
the  increased  current.  In  the  outer  regions  of  the  halo,  heat¬ 
ing  leads  to  an  expansion  of  the  gas;  towards  the  center  of  the 
halo  heat  conduction  must  become  important  as  the  hot  gas  becomes 
part  of  the  arc-like  central  column  of  the  streamer. 


A  detailed  analysis  must  include  all  these  interacting 
nonlinear  processes  in  a  consistent  manner.  We  do  not  yet  have 
such  a  description,  but  we  have  made  substantial  progress  toward 
it.  This  year  we  have  included  in  our  model  gas  heating  ahead 
of  the  streamers,  near  the  edge  of  the  halo,  and  within  the  halo, 
making  it  possible  to  estimate  the  relative  importance  of  these 
different  effects  in  the  different  regimes  where  the  discharge 
can  be  operated.  Moreover,  we  now  have  what  we  believe  to  be  a 
more  precise  description  of  the  radial  variation  of  vapor  den¬ 
sity,  electric  field  and  conductivity  within  the  halo,  a  region 
which  we  had  previously  been  able  to  describe  only  phenomenologi - 
cally.  By  assuming  that  the  field  within  the  halo  is  damped  at 
the  breakdown  threshold,  and  combining  this  with  an  analysis  of 
the  vapor  expansion  due  to  heating,  we  have  been  able  to  derive 
scaling  laws  which  put  our  earlier  model  on  a  firm  foundation 
and  suggest  a  way  in  which  it  should  be  possible  m  future  work 
to  construct  from  first  principles  a  more  complete  description. 

In  what  follows,  our  earlier  calculations  are  reviewed  and  then 
the  more  recent  work  is  described.  Since  the  earlier  work  was 
covered  in  last  years 's  final  report  *  ,  we  are  including  here 

only  a  brief  summary  of  the  essential  points. 

2.  STREAMER  CURRENT  AND  CORE  RADIUS 

Since  the  electrical  conductivity  is  high  within  the 
streamer,  the  streamer  current  is  determined  by  the  resistivity 
of  the  background  discharge.  The  streamer  is  essentially  a 
resistively  ballasted  arc,  with  the  cold,  e-beam  ionized  gas 
ahead  of  the  streamer  imposing  the  limiting  resistance.  Since 
the  e-beam  generated  conductivity  is  uniform,  and  the  current  is 
quasi-static  (i.e.,  divergenceless)  one  has, 

V  (a  E )  =  o  V.  E  =  -  •  V‘V  -  0 

*  o—  o  •  -  o 


So  the  potential  obeys  Laplace's  equation  and  hence  is  conveniently 
described  in  spherical  coordinates  in  terms  of  Legendre  polynomi¬ 
al  1  s  .  Assuming  that  the  streamer  halo  is  a  conducting  sphere  with 
uniform  surface  potential  Vj ,  that  the  E-field  far  from  the 
streamer  is  uniform  and  of  magnitude  E0,  and  that  the  potential 
at  the  base  of  the  streamer  (f*  =  -,  r  =  c)  is  zero,  one  obtains 
a  reasonably  accurate  description  in  terms  of  the  first  few 
Legendre  polynomials: 


V  =  E  (r 
o 
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5.  Electron-Beam  Sustainer  Discharge  Streamers  and  Arcing,"  Final 
Technical  Report  (Nov.  1977). 
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where,  again,  c  is  length  of  the  streamer  column,  and  a  is  the 
radius  of  the  halo  or  cap  sphere.  Since  the  conductivity  of  the 
streamer  is  much  higher  than  that  of  the  background  gas,  the  po¬ 
tential  Vj  of  the  halo  surface  must  be  close  to  zero.  So  at 
least  as  a  first  approx imat ion  (which  is  probably  as  accurate  as 
our  assumption  that  the  halo  is  at  uniform  potential)  one  can  set 
V}  0.  Then  differentiating  Eq .  (2)  to  find  the  strength  of 

the  radial  field  at  the  edge  of  the  halo  gives. 


E 


r  =  a 


E  [(1  +  -  ♦  — )  +  3  cos  <•) 
o  a  c 


(  J) 


From  this  equation  one  can  determine  the  streamer  current. 
Setting  j  *0E  and  integrating  over  the  surface  of  the  cap 
sphere  gives  the  total  current, 


I 


/ 
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(4) 


where  the  second  term  i n  Eq .  (3)  has  been  integrated  over  the 
upper  hemisphere  and  a/c  has  been  neglected  because  it  is  much 
smaller  than  c/a. 

Our  original  model  also  used  Laplace's  equation  to  esti¬ 
mate  the  radius  of  the  top  of  the  streamer  column,  i.e.,  of  the 
hot  core  inside  the  halo.  Reasoning  that  while  the  conductivity 
is  not  exactly  constant  within  the  halo  (it  probably  does  not 
change  too  abruptly),  one  can  estimate  b,  the  core  radius,  by  tak¬ 
ing  V  =  0  at  r  -  b  and  expanding  in  Legendre  polynomials  as  above. 
By  then  assuming  that  E  =  E0  at  r  =  a  one  obtains  an  estimate  of 
the  ratio  a/b. 


3.  TEMPERATURE  OF  THE  STREAMER  COLUMN 

In  the  main  column  of  the  streamer,  the  electron  tempera¬ 
ture  and  density  are  in  equilibrium  with  the  neutral  gas,  so  the 
degree  of  ionization  obeys  Saha's  equation  at  the  temperature  of 
the  neutral  gas.  Tt>e  electric  field  is  low,  but  the  current  den¬ 
sity  is  high,  so  the  resistive  cower  input  is  sufficient  to  keen 
the  column  hot.  In  gas  at  atmospheric  pressure  with  a  tempera¬ 
ture  below  3  x  lO-*  °K,  heat  loss  is  mainly  by  conduction. 
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Without  knowing  the  details  of  the  column 
solve  the  heat  flow  equation  exactly,  but 
proceeded  by  averaging  it  over  the  column 


structure  one  cannot 
following  Raiser^’*  we 
area.  This  gives. 


o 


A  K  T 


(6) 


where  b  is  the  column  radius.  T  is  the  mean  temperature  of  the 
gas  within  the  column,  <  is  the  thermal  conductivity,  and  A  is 
a  numerical  factor  characteristic  of  the  gas.  For  high  tempera¬ 
ture  air,  Kaiser's  value  of  A  is  15,  which  we  used  for  nitrogen 
in  our  model.  The  mean  thermal  conductivity. 


T 

ic  (T'  )  dT  ' 

o 

(  7 ) 

is  derived  from  the  published  conductivity  for  air  at  high 
temperatures.  Over  the  range  of  interest,  1500°K  <  T  <  1500°K, 
•  can  be  approximated  by  the  analytic  expression, 


K  T1*893 
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3.145  x  10“10  T1*893  W/cm2°K 


l  7 ) 


Under  pressure  balance,  the  gas  density  is  inversely  pro¬ 
portional  to  temperature,  and  the  electrical  conductivity  scales 
directly  with  electron  density  and  inversely  with  gas  density,  so 


T  ,  /ne<T>  *  neo  V 

T  1  (n  IT  )  +  n  /  ° 

o  \  e  o  eo/ 
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where  0,  ro  and  neo  are  the  conductivity,  temperature,  and  e-beam 
generated  secondary  electron  density  of  the  background  discharge. 
The  electron  density  in  the  streamer  column,  ne(T)  is  described 
by  the  Saha  equation. 


where 
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and  L0  is  Loschmidt's  number,  Pc  is  the  pressure  in  atmospheres, 
ui  and  uc  are  the  statistical  weights  of  the  ion  and  atom.  6  is 
the  degree  of  dissociation  of  a  molecular  gas  (6=2  for  nitrogen 
in  the  regime  of  interest!  me ,  k,  and  h  are  the  electron  mass, 


6.  Raiser,  Y.P.,  Sov.  Phys.  JETP,  31,  1148  (1970). 

7.  NACA,  TN  4150  (1958) . 
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Boltzmann  constant,  and  Planck  constant,  respectively,  and  1  is 
the  ionization  potent  1  a  1 ( a '  . 

Combining  Eqs.  (4)  through  (8)  one  obtains, 


Together  with  Eq .  (9)  for  ne(T)  this  result  gives  the  tempera¬ 

ture  T  of  a  thermally  stable  streamer  column.  The  result  does 
depend  upon  b,  but  only  very  weakly.  For  a  typical  discharge 
in  nitrogen  (c  1cm,  b  -  1mm,  T0  =  300°K,  'Q  =  2.6  x  10“4 
mho/cm,  E0  =  5000  V/cm) .  These  equations  predict  that  T  =  6200°K. 
While  some  of  the  equations  used  in  this  calculation  are  merely 
rough  estimates,  the  variation  of  electron  density  with  tempera¬ 
ture  is  so  strong  that  temperature  predicted  by  this  calculation 
should  be  reasonably  accurate. 


4.  REQUIRED  ENERGY  INPUT 


We  have  hypothesized  the  streamer  growth  mechanism  to  be 
resistive  heatinq  by  the  concentrated  current  flowing  into  the 
streamer,  its  temperature  rises  toward  that  of  the  streamer 
column,  making  it  part  of  the  streamer,  and  causing  the  current 
concent  ration  to  extend  farther  into  the  background  discharge. 

To  analyze  this  process,  one  must  know  the  energy  required  to 
heat  the  gas  to  the  column  temperature  predicted  by  the  calcula¬ 
tion  described  in  the  preceding  section. 

In  the  subsonic  regime,  the  energy  required  to  raise  one 
gram  of  gas  at  constant  pressure  from  temperature  T0  to  T  is 
equal  to  the  enthalpy  difference  between  these  temperatures. 
Writing  enthalpy  as  q(T).  the  energy  is  AQ  =  q(T)  -  q(To>  joules 
per  gram.  The  energy  required  to  raise  unit  volume  of  gas  over 
the  same  temperature  range  at  constant  pressure  will  then  be, 


Zeldovich,  Y.B.,  and  Raizer,  Y.P.,  Physics  of  Shock  Waves  and 
High  Temperature  Hydrodynamic  Phenomena,  A.P.  New  York  (1966) . 
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where  .  =  ,  ( T * )  is  the  density  of  the  gas.  Since  t  -  ,0  T0/T, 
where  is  the  density  at  T  ,  this  reduces  to 


T 

o 


dT’ 

T' 


Physically,  as  a  constant  volume  of  gas  is  heated  under  constant 
pressure,  the  gas  expands  and  some  leaves  the  volume,  removing 
energy,  but  leaving  less  gas  in  the  volume  to  be  heated  further. 
The  above  expression  for  •  ,  the  "constant  volume  constant  pres¬ 
sure  enthalpy"  includes  these  effects. 

For  molecular  gases.  (T)  can  be  quite  comolicated,  but 
over  a  limited  ranae  of  temperature  it  can  usually  be  approxima¬ 
ted  by  a  simple  power  law.  For  nitroqen  in  the  region  of  inter¬ 
est  to  us,  data  available  in  the  thermodynamic  tables<9)  1S  ade¬ 
quately  fitted  by  the  simple  formula. 


e  ceo 


.6  1 


where  tQ  -  1.0  J/crr  .  Tj  =  3500oK,  Po  is  the  oressure  in  atmos¬ 
pheres,  and  t  is  the  energy  required  to  heat  1  cm -3  of  gas  from 
30  0 3  K  t  o  T .  For  t :  •  streamei  Lumn  temper  a  tun  :alcuiated  il  ve , 
6200=’K,  this  gives  a  needed  input  energv  of  -1.4  J/cm^. 


b.  HEATING  AT  THE  HALO  EDGE 


In  our  original  model,  we  calculated  the  resistive  heating 
within  the  halo  under  the  assumption  that  the  conductivity  varied 
radially  as 


1  ) 


This  form  w  is  simply  postulated,  but  the  heatino  rate  it  qivcs 
should  be  a  reasonable  estimate  if  the  conductivity  does  increase 
rapidly  as  one  moves  into  the  halo.  (In  support  of  this  we  con¬ 
sidered  alternate  forms  and  showed  that  the  result  was  practi¬ 
cally  unchanged.)  Using  this  expression  one  can  integrate  j 2 / , 
over  the  cap  sphere  from  the  edge  of  the  core  to  the  edge  of  the 
halo  to  obtain. 
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Since  the  ratio  of  conductivities  in  the  core  at  6200°K  and 
in  the  background  gas  under  e-beam  ionization  is  appnmately 
10^  and  a/b  x  3,  the  value  of  k  which  matches  the  boundary  con¬ 
dition  on  the  conductivity  is  k  =  3.45.  Thus,  the  second  term 
in  the  above  equation  is  completely  negligible  and  the  mean 
power  density  within  the  halo  is 

3 1  2  . .  ,  3 

q  =  - x - 2 -  W/cm 

16  r  a  k 

o 


From  this  power  density  one  could  compute  the  time  required 
to  heat  the  halo  to  the  temperature  of  the  core  of  the  streamer. 
However,  a  heating  time  calculated  in  this  way  could  be  too  long 
because  only  the  center  need  be  heated  for  the  core  to  grow.  In 
a  calculation  given  in  detail  in  the  previous  annual  report  we 
deduced  a  correction  factor  of  '  =  0.3,  which  accounts  foi 
radial  temperature  variation  consistent  with  the  above  stated 
radial  dependence  of  the  conductivity. 

From  these  results  one  can  deduce  a  character ist ic  velocity 
for  propagation  of  the  streamer.  It  is  the  radius  of  th<  sphere 
a,  divided  by  the  time  required  to  heat  the  gas  in  the  halo  in¬ 
terior  to  the  core  temperature,  '  -  '*  /q.  Thus  we  have, 
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where  we  have  used  Eq. 

(4) 

to 

express 

the 

current  I  in  terms  of 

the  other  variables. 

The  model  predicts  that  c/a  should  remain  roughly  constant 
as  the  streamer  -rows  and  this  is  confirmed  by  the  framing  camera 
photoqraphs.  Putting  the  result  in  a  form  more  appropriate  for 
such  a  scaling,  and  using  the  above  derived  values  of  ,  k,  and 
,  we  have 


where 
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Hence  the  model  predicts  that  streamei ;  should  qrow  exponentially. 
This  prediction  is  in  qood  agreement  with  our  observations.  The 
predicted  growth  rate,  .  is  primarily  dependent  upon  the  speci¬ 
fic  power  loading,  'o^'o^/Pc'-  This  is  also  in  good  agreement  with 
the  data  over  a  considerable  range  of  parameters.  As  can  be  seen 
from  Figure  5,  the  line  of  constant  growth  rate  is  nearly  a  line 
of  constant  power  density  in  the  stronq  field  regime.  Since  the 
assumptions  underlying  our  theory  (e.g.,  that  the  halo  around  the 
streamer  tip  extends  well  beyond  the  radius  of  the  streamer 
column)  are  ippropriate  to  the  strong  field  operating  mode,  one 
expects  agreement  with  experiment  there.  In  the  high  current, 
low  field  regime,  the  resistive  heating  of  the  gas  beyond  the 
halo  may  be  important  and  that  heating  was  not  included  in  the 
riginal  model.  We  have  done  some  analysis  of  heating  ahead  of 
the  halo,  as  described  in  the  next  section,  but  in  the  regime 
of  principal  interest  for  laser  appl icat ions ,  our  original  model 
is  in  good  agreement  with  the  data. 

t.  HEATING  AHEAD  OF  THE  STREAMER 

The  assumed  variation  of  conductivity  within  the  halo  re¬ 
stricted  our  original  model  to  effects  of  heating  just  inside  the 
halo  edge.  Clearly,  a  more  detailed  theory  should  include  a  cal¬ 
culation  of  the  power  input  throughout  the  vicinity  of  the  grow¬ 
ing  streamer.  We  have  now  made  considerable  progress  toward 
the  development  of  a  more  complete  descr lpt ion . 

There  is  some  heating  beyond  the  halo,  which  can  be  esti¬ 
mated  ft  tn  •  he  11  ilysi  ilri  i  ly  lone.  fust  .t  .  :•  t  n<  •  :  • 

electric  field  is  Ec,  the  critical  field  and  the  conductivity  is 
-0.  Therefore,  the  power  input  is  o F.c2,  which  will  heat  the 
gas  to  the  conductive  temperature  in  a  time  :  =  >  -0EC2.  Since 

*  he  scale  length  of  this  field  is  a,  the  halo  radius,  such  heat¬ 
ing  alone  would  give  a  velocity  v  =  a/*.  We  have  not  yet  dir¬ 
ectly  estimated  a,  but  it  is  implied  by  Eq.  (3)  for  the  E-field 
beyond  the  halo.  Setting  •>  =  90r,  and  F  Ec  we  have 
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This  is  approximate,  because  the  ••  dependence  has  been  omitted, 
but  is  consistent  with  our  original  approximation  of  the  halo  as 


a  sphere  (which  was  used  to  derive  Eq. 
the  propagation  speed 
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From  this  we  have 
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This  is  primarily  dependent  upon  ’o  E0  -  ]0>  the  current  density 
within  the  background  discharge.  This  is  an  extreme  conclusion, 
since  it  ignores  all  heating  within  the  halo,  but  we  note  that  at 
low  electric  fields  where  heating  beyond  the  halo  becomes  more 
important,  the  data  in  Figure  5  are  not  inconsistent  with  a  con¬ 
stant  growth  rate  line  being  a  line  of  constant  current  density 
in  the  weak  field,  high  current  regime. 

7.  THE  STRUCTURE  OF  THE  HALO 


We  believe  we  have  now  found  a  way  to  describe  the  whole 
region  around  the  streamer  tip  in  more  detail.  Our  new  approach 
starts  with  the  assumption  that  inside  the  halo,  but  outside  the 
core  cylinder,  the  electric  field  strength  is  clamped  at  the  break¬ 
down  threshold  level.  We  reason  that  the  current  is  limited  by  the 
background  discharge  and  hence  cannot  increase,  so  the  E-field  must 
be  consistent  with  this  current.  If  the  E-field  dropped  below  E., 
the  conductivity  would  drop  causing  the  field  strength  to  increase 
again.  But  if  the  field  rose  much  above  E  ,  the  conductivity  would 
be  rapidly  increased  by  electron  multiplication,  which  would  reduce 
the  electric  field  toward  Ec .  We  thus  conclude  that  in  the  halo 
the  field  strength  cannot  be  very  different  from  the  breakdown 
field. 


However,  this  does  not  mean  that  the  field  strength  is  con¬ 
stant  in  the  halo,  because  t he  breakdown  field  depends  upon  the 
gas  density,  which  decreases  as  the  gas  is  heated.  More 
precisely,  _ 

c.  co 

E  =  - 

n 

o 

where  by  Eco  is  meant  the  breakdown  field  is  the  background 
discharge. 


Within  the  halo  one  expects  the  current  and  the  E  field  to 
be  nearly  radial,  since  both  are  diverging  from  the  tip  of  the 
streamer  column.  This  implies  that  the  magnitude  of  the  current 
density  scales  as 
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Hence  the  power  density  must  scale  as 

nE 

r.  t.  co 

P  -  j.t  *  — , 

n  i  *" 

o 

Matching  this  to  the  boundary  condition  at  the  halo  edge  gives 
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For  nitrogen,  it  was  demonstrated  earlier  that 
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Differentiating  this,  we  have 
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Combining  these  results  gives 
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Now  make  the  substitution  dn/dt  v  dn  di  t  ft  air. 
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This  has  the  solution,  for  n  =  n  at  r  =  a, 
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Well  inside  the  halo  this  reduces  to 


whe  re 
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This  result  provides  the  scaling  needed  to  develop  a  more 
comprehensive  theory.  That  has  not  yet  been  done,  but  we  note 
that  the  new  result  does  provide  a  more  solid  foundation  for  our 
earlier  model.  From  the  above  results  we  have  that  power  density 
well  inside  the  halo  varies  radially  as 
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Since  this  increases  with  decreasing  r,  the  power  density  is 
highest  near  the  center  of  the  halo,  but  since  the  increase  is 
less  rapid  than  l/r2,  most  of  the  total  power  input  to  the  halo 
is  near  the  outer  edge  as  we  assumed.  However,  we  now  see  that 
heating  near  the  center  is  not  negligible  and  must  be  included 
in  a  more  complete  description.  (However,  we  have  not  yet  in¬ 
cluded  heat  conduction,  which  will  further  reduce  the  power  in¬ 
put  near  the  core.) 

Clearly,  this  last  calculation  is  just  the  first  step  in 
the  development  of  a  second  generation  theory,  but  we  think  it 
is  important  because  it  shows  how  to  proceed  with  that  develop¬ 
ment  . 


MAGNETIC  DEFLECTION  OF  THE  STREAMER  COLUMN 


In  the  course  of  this  year's  experiments,  a  magnetic  field 
was  imposed  transverse  to  the  discharge  to  examine  the  possi¬ 
bility  that  this  would  push  the  streamers  sideways  through  the 
cold  background  gas.  No  such  effect  was  seen,  although,  as 
discussed  above,  the  B-field  did  seem  to  affect  the  streamer 
location.  The  strongest  field  compatible  with  our  e-bcam  was 
employed  ~200  gauss,  which  gives  a  100  kV  electron  a  gyro  radius 
of  5  cm,  comparable  to  the  dimensions  of  our  e-beam.  But  a  simple 
calculation  shows  that  even  this  was  marginal  for  the  intended 
purpose.  The  current  density  in  the  background  gas  in  our  dis¬ 
charge  is  typically  1.5  amperes/cm^.  Since  a  wel 1 -devel oped 
streamer  draws  the  current  from  about  1  cm^  of  discharge,  we 
know  that  a  typical  streamer  current  is  around  1.5  amperes.  A 
magnetic  field  of  200  gauss  would  thus  impose  upon  the  streamer 
column  a  force  of  about  30  dynes/cm.  A  streamer  column  is  a 
fraction  of  a  millimeter  in  diameter,  and  the  gas  density  is 
reduced  by  at  least  an  order  of  magnitude  from  atmospheric  den¬ 
sity,  so  a  reasonable  estimate  of  the  mass  of  the  gas  in  the 
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column  is  10  g/cm.  If  there  were  no  resistance,  the  maqnetic 
force  would  move  this  mass  a  millimeter  or  two  in  the  few  tens 
of  microseconds  available. 

We  would  have  observed  such  motion,  but  since  there  is  re¬ 
sistance  from  the  gas  surrounding  the  column  it  is  not  surpris¬ 
ing  that  the  streamers  did  not  move  this  far.  Streamer  motion 
sideways  through  the  discharge  is  a  complicated  process,  and 
there  appears  to  be  no  simple  way  to  predict  exactly  how  much 
B-field  is  required.  We  thought  there  was  a  possibility  of  ob¬ 
serving  some  small  deflection  m  our  apparatus,  which  would  have 
given  us  a  data  point  from  which  we  could  extrapolate,  but  this 
did  not  occur.  Since  stronqer  fields  are  not  usable  in  our 
apparatus,  this  precludes  further  investigation  of  the  process 
there . 

However,  there  remains  a  distinct  possibility  that  an 
external  B- field  would  affect  streamers  in  other  discharge 
machines.  One  could  apply  a  stronger  B-field  parallel  to  the 
e-be.un.  Indeed  this  is  done  in  some  larqe  machines  to  minimize 
e-beam  spreading.  Since  the  sustainer  discharge  is  normally 
parallel  to  the  e-beam,  such  a  field  would  not  initially  exert 
a  force  on  streamers.  But  we  know  that  in  flowing  gas  systems, 
streamers  are  swept  downstream,  where  they  lead  to  downstream 
arcinq.  A  streamer  extending  downstream  would  be  transverse  to 
the  B-field,  so  there  could  be  an  effect.  Moreover,  the  times 
involved  are  lonqer  than  in  our  machine,  and  the  currents  are 
larger,  so  our  observations  in  no  way  preclude  the  possiblity  of 
large  deflections  and  resultant  cooling  of  such  streamers. 


SECTION  IV 


CONCLUSIONS 


In  summary,  extensive  new  framing  camera  data  on  streamer 
structure  and  development  have  well  confirmed  our  model  of  a 
streamer  as  composed  of  core  and  halo  surrounded  by  background 
discharge.  Streamer  growth  is  consistently  seen  to  be  exponen¬ 
tial,  with  the  ratio  of  length  to  halo  radius  remaining  constant 
during  growth  for  a  wide  range  of  discharge  conditions.  Our  orig¬ 
inal  theory  predicts  a  growth  rate  primairly  dependent  upon  power 
loading,  and  we  find  this  to  be  true  in  the  high  field  regime 
where  laser  discharges  are  operated,  but  not  in  the  weak  field, 
high  current  regime.  The  model  has  been  further  developed  by 
calculation  of  heating  ahead  of  the  streamer  and  in  the  interior 
of  the  halo.  Scaling  laws  for  the  density,  temperature,  and  elec¬ 
tric  field  distributions  within  the  halo  have  been  derived,  and 
these  provide  a  basis  for  the  construction  of  a  more  complete 
theory  in  future  work. 

This  year's  work  has  established  the  need  for  a  detailed 
mapping  of  the  gas  density  and  resistive  power  distributions 
around  the  tip  of  a  streamer.  We  believe  this  can  be  done  theo¬ 
retically,  and  we  have  also  proposed  usinq  shadowgraph  techniques 
to  obtain  data  on  the  density  variation. 

The  imposition  of  a  transverse  magnetic  field  was  not  seen 
to  affect  the  streamer  growth  rate  or  to  move  the  streamer  later¬ 
ally,  but  practical  limitations  of  our  particular  discharge  pre¬ 
vented  a  complete  examination  of  that  question. 

This  year’s  work  provided  the  first  framing  camera  records 
of  spontaneous  streamer  development  from  a  smooth  cathode  surface. 
Considerable  scatter  is  seen  in  the  streamer  formation  times,  in¬ 
dicating  that  more  detailed  study  of  the  initiation  process  in  the 
cathode  layer  would  be  valuable.  This  will  also  be  pursued  in 
future  work . 
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